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Abstract 

The Kardar-Parisi-Zhang (KPZ) equation is accepted as a generic descrip- 
tion of interfacial growth. In several recent studies, however, values of the 
roughness exponent a have been reported that are significantly less than that 
associated with the KPZ equation. A feature common to these studies is the 
presence of holes (bubbles and overhangs) in the bulk and an interface that 
is smeared out. We study a model of this type in which the density of the 
bulk and sharpness of the interface can be adjusted by a single parameter. 
Through theoretical considerations and the study of a simplified model we de- 
termine that the presence of holes does not affect the asymptotic KPZ scaling. 
Moreover, based on our numerics, we propose a simple form for the crossover 
to the KPZ regime. 

PACS numbers: 05.40.-a; 05.70.Np; 68.35.Rh 
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Since its inception over fifteen years ago, tlie Kardar-Parisi-Zliang (KPZ) equation |T| 
lias proven itself as a generic stochastic description of a roughening interface. Such inter- 
faces arise in a vast range of physical situations ranging from colloidal aggregation through 
bacterial growth to forest fires 0. Typically one describes such systems in terms of a d- 
dimensional plane (substrate) above each point r of which one associates an interfacial height 
function /i(r). The time evolution of the surface described by the height function is given 
by the rf+l dimensional KPZ equation thus: 

^h(v, t)=v^ + uV^h + ^(V/i)2 + 77(r, t) . (1) 

Here ^^(r, t) is a Gaussian white noise of zero mean and correlator (77(r, t)ri{r', t')) = F^'^dr — 

r'm-t')- 

Only for a small number of microscopic growth models — such as restricted solid-on-solid 
(RSOS) models — has (|l|) been obtained from first principles p. Nevertheless, the ap- 
propriateness of (|I]) as a coarse-grained description is easily justified using phenomono logical 
arguments as follows. Firstly, the constant Vq is simply the mean rate at which a fiat in- 
terface would proceed in the absence of noise; the Laplacian represents a tendency of the 
interface to smoothen through surface tension and the non-linear term expresses the fact 
that one generally expects growth to occur normal to an interface. The noise term models 
fluctuations superposed on deterministic growth rules. 

Kinetic roughening is conveniently studied through examination of the interfacial width 
Wit): 

W\t) = {Ah{ty) ; Ah{ty = hi^ - Hff . (2) 

The angle brackets represent an average over different realizations of the noise (an ensemble 
average) whereas the overbar denotes a spatial average in a given realization at a given 
time t. It has long been known that the behavior can be summarized as a dynamic 
scaling relation W{t) ~ t^f{t/L^) where /(M)=const for small u, giving j3 as the early time 
growth exponent; for large times t ^ L^, the width saturates as 1^ ~ L° and f{u) ~ u"^ 
implying that a=l3z. For the 1+1 dimensional KPZ universality class the exponents are 
known exactly (via a range of methods, see 0) as a=|, /?=| and z=j. 

The robustness of the KPZ description is due to a very important property: the deriva- 
tives of h that appear in (|T]) are the only ones whose contribution does not diminish under a 
rescaling transformation. The consequence of this is that any interface whose rate of growth 
depends only its local shape will be adequately described by (|T]) on sufficiently large length 
and time scales. Known ways of changing the asymptotic scaling are to alter the nature 
of the noise — e.g. to power-law distributed or correlated noise — or to introduce nonlocal 
dynamics. By nonlocal dynamics it is meant that the velocity of the interface at a point 
depends on quantities other than the gradient of the height function there. 

One origin of nonlocal contributions to surface growth is a bulk structure that is not 
compact. To understand this, consider first the case where the structure is compact, i.e. the 
surface grows by absorption of particles in such a way that no holes (bubbles or overhangs) 
are created. RSOS models are examples of such a system. In this instance it is clear 
that the interface is sharp on the microscopic scale and the entire system can be completely 



2 



described by the height function h{r). Furthermore, the interfacial dynamics are entirely 
local and equation ([l|) should apply. Indeed it is found that RSOS models exhibit KPZ 
scaling even at small system sizes 0. 

An example of the contrasting situation of a noncompact bulk structure is the Fisher wave 
1^. To realize a Fisher wave microscopically, one must include particle removal processes in 
a growth model. As a result: (i) the bulk contains holes; (ii) the interface is smeared out: 
by this we mean that there is a finite interfacial region in which a coarse-grained density 
field decays to zero due to the presence of holes; and (iii) the interfacial motion could, 
in principle, be affected by fluctuations within the density field away from the interface 
(nonlocal dynamics). Consequently, the height function h{r) does not uniquely specify 
a configuration of the density field and thus it is possible that its evolution may not be 
adequately described by the KPZ equation (|1|) which contains only h{r). Numerical evidence 
for this scenario has recently been presented in the context of Fisher waves: most notably 
through a roughness exponent of a^OA for d=l [^]. 

Even when the interfacial dynamics are local, numerical studies have suggested that a 
is reduced when the bulk structure contains holes: for the ballistic deposition (BD) model 
values of a=0.42(3) and 0.47(1) Q were reported in early work on 1+1 dimensional 
BD, whereas more recent very large-scale simulations ||^ yielded an estimate of a=0.45. 
Although BD is generally accepted as being a realization of KPZ growth, the possibility 
that a noncompact bulk structure implies a new non-KPZ university class has been raised 
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In this work we introduce a model of interfacial growth in which a parameter k allows 
us to vary the density of the bulk, hence the sharpness of the interface and in turn any 
nonlocality in the interfacial dynamics. When k=0, the bulk is essentially compact and the 
interface nearly sharp — i.e. the density decays rapidly to zero across the interface. Here, we 
expect, and indeed observe, KPZ interfacial scaling. For nonzero k (a more smeared-out 
interface) our preliminary studies give an estimate of alpha close to 0.4, consistent with [Q. 
Thus the model apparently gives two different scaling behaviours depending on whether the 
interface is smeared out or not. However, we argue, through consideration of how the KPZ 
equation might be modified to incorporate the effect of holes in the interfacial region, that 
after sufficient coarse-graining the KPZ equation is in fact the correct description for all 
values of k. In order to demonstrate this quantatively we study a simplified model that 
allows us to push the numerics to larger systems. We find that our data can be well fitted 
by a crossover scaling form to the KPZ scaling. The considerations that lead us to conclude 
that the KPZ equation is the correct description are generic and not model specific. Thus we 
conclude that values of a roughness exponent close to 0.4, as reported elsewhere for different 
models are very likely the result of not having reached the asymptotic scaling regime. 

We now define our interfacial growth model, which we shall call the wetting model. 
It comprises a c?- dimensional substrate at z=0 where z is the growth direction. Periodic 
boundary conditions (with periodicity L) are imposed in the directions perpendicular to z. 
In order to sustain growth, the substrate is kept fully occupied (wet) at all times. Wetting 
events, where an occupied site causes its neighbor above (in the +z direction) to become 
occupied, occur at a rate r. Additionally, a site at z>0 can dry out (become empty) with 
rate k. Particles (occupied sites) can move to a neighboring site with the same z co-ordinate, 
each at rate D/{2d) and subject to the condition that the receiving site is empty. We define 
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the height h of the interface between the wet and dry regions as the z co-ordinate of the 
uppermost occupied site above substrate position r [|ll|]. We consider henceforth only the 
case d=l (see Fig. |l|) and thus replace r with x. 

Before concentrating on the interfacial scaling behavior, we note some other interesting 
properties of the wetting model. Firstly, if the ratio r/k of the wetting rate to the drying-out 
rate is too small, there will be no interfacial growth: as r/k is increased past some critical 
value (dependent on D) there is a wetting transition [|l^ to a regime where the wet region 
invades the dry part of the lattice. This transition is related to the directed percolation 



universality class [|13|. Under these conditions one has a moving interface and interfacial 
growth. The limit D=0 reduces the model to a set of L non-interacting Id stochastic 
processes, namely asymmetric contact processes. The opposite limit D ^ oo, L ^ oo 
renders the dynamics equivalent to a deterministic version of the same Id process. Indeed, 
this case constitutes a microscopic realization of deterministic growth processes of the type 
studied in A fuller discussion of all the regimes will be presented elsewhere for the 
present study of interfacial scaling, we use intermediate values of D and r. 

We first check the growth of the interface in the case k=0 which is the situation where no 
holes may be spontaneously created in the bulk. Thus the bulk density is one, and the only 
holes present are in a small zone (numbering a few lattice sites) near the interface and which 
have their origin in the sideways diffusion process. We consider the interface therefore 
effectively sharp and thus expect to find KPZ behavior. This is tested by performing a 
scaling collapse of the width {W{t)/L'^ against t/L^) which is shown in Fig. |^(i) for a range 
of system sizes and representative values of D,r. Good data collapse was achieved using 
exponents consistent with KPZ scaling: 2;=!. 50(5), a=0.49(2) where the figure in brackets 
indicates the change in the last quoted digit required before the collapse becomes noticeably 
worse. 

The second plot in Fig. ^ shows the corresponding collapse for the case where k=l, i.e. 
wet sites can dry out leaving holes in the bulk. We found we had to change the scaling 
exponents to z=1.35(5), a=0.41(2) in order to obtain collapse of the saturation width. Note 
that in comparison to Fig. 0(i), the collapse is far less convincing. Therefore, although we 
appear at first sight to have obtained a value of the roughness exponent a consistent with 
MM, it is not clear that we were able to probe the true asymptotic scaling regime before 
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FIG. 2. The width function W{t) in the wetting model with the drying-out process (i) inactive 
and (ii) active. Data from different system sizes L were made to collapse (at least at late times) 
by scaling W with L° and t with . The values of a, z used in each case are shown on the plots. 



simulation run times rendered increases in L impractical. 

To assess whether the true scaling regime had been accessed, i.e. whether one truly has 
asymptotic scaling different to KPZ in figure H(ii), we consider how a smeared out interface 
might be described through a modifed version of (|I]). In order to incorporate a generic 
coupling of the interface to holes in the bulk we introduce an additional term into (P as 
follows 

t) = Vo + lyVh + ^{^hf - kg{x) + r] (3) 

in which we have included an explicit dependence on k, so that when k=0 this equation 
reduces to (|l|) as required. The function g{x) is intended to capture any nonlocal dynamics 
present. In the case of the wetting model, —kg{x) is the rate at which h{x) decreases when 
a particle at the interface disappears, implying that g{x) is the size of the gap behind the 
interface, defined as the distance between the uppermost and second uppermost particle at 
position X. We have checked explicitly that decreases in h{x) follow closely the profile of the 
gap size g{x) and so we believe that describes the wetting model adequately, at least on 
a coarse-grained level. 

We now proceed to show that the presence of the extra term in @ does not affect the 
long- wavelength properties of the interface. To this end we studied the statistics of g{x) 
in the wetting model. We found first that the distribution of gap sizes became stationary 
rather quickly, at times when the interface was still roughening. In Fig. ^(i) we plot the 
stationary gap size distribution for three different system sizes. We note that the distribution 
is system-size independent and decays exponentially. Thus the finite length scale associated 
with the gap remains constant as the substrate length is increased, and is irrelevant once L 
is sufficiently large. Of course, at smaller L, g{x) will play some role. 
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FIG. 3. (i) The gap-size g probability distribution for the wetting model with L=100, 200, 400. 
(ii) The gap-size correlation function in the x-direction for the wetting model with L=100, 200, 400. 
In both cases D=5, r=1.5 and k=l. 



Now we argue that unless g{x) exhibits scale-invariant correlations in the x direction, it 
will be rescaled into the noise term already present in (|l|). Specifically, g can be replaced 
with go + fj where go is a constant which can be absorbed into vo and is a noise term 
with irrelevant short-range correlations that can be absorbed into the noise term 77 in (P. 
To exclude the possibility of long-range correlations, we plot the correlation function c{6) = 
{g{0)g{5)) /{g{OY) - 1 in Fig. |(ii). It is clear that for L=100, 200, 400 the correlation length 
retains the same (finite) value. The small anticorrelation for large S vanishes with increased 
L and is a consequence of finite system sizes. 

Thus, our results strongly suggest that the term g{x) modifying the KPZ equation from 
(H) to d^) is, for the wetting model, purely cosmetic. In short, we contest that if L were 
increased beyond the values used in Fig. ^ a crossover towards KPZ scaling would be ob- 
served. However, the system sizes that would be required to demonstrate quantitatively the 
crossover in the wetting model are unfeasible. 

In order to examine the nature of such a crossover we study a simplified model which 
retains the important feature of the wetting model, namely a parameter that allows us to 
go from a sharp to a smeared-out interface. The model is constructed by replacing the 
two-dimensional density field of the wetting model with an interface h{x) coupled to a gap 
of size g{x). Thus this new model, which we will refer to as the ballistic deposition and 
desorption (BDD) model, should also be governed by (^. 

The dynamics of the BDD model are defined as follows. At each time step, one chooses 
a substrate position x at random then performs either a particle deposition move (at unit 
rate) or a desorption move (at rate k). In the former case, a particle is "dropped" vertically 
downwards in column x until it comes to rest at a site whose nearest-neighbor is occupied. 
This increases the height h{x) by an amount which defines the gap size g{x). A desorption 
move is implemented by decreasing h{x) by g{x) and replacing g{x) with one of the other 
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gap sizes in the system, chosen at random (maintaining a self-consistency in the gap-size 
distribution). The desorption move serves to smear out the interface, and thus the rate k 
plays the same role as k in the wetting model (although no numerical equivalence between 
the two should be assumed). 

With the desorption rate k set to zero, one recovers the ballistic deposition model which 
is a widely accepted realization of (P [|l6l and thus the interface should exhibit KPZ scaling 
behavior. With nonzero k the height h can decrease by a random variable g. We have found 
that the statistics of g in the BDD model are very similar to those shown in Fig. ^ for the 
wetting model: both the gap distribution and correlation functions are stationary and have 
no dependence on system size; the gap size distribution has an exponential tail and the 
correlation length in the x direction is almost zero. Thus we believe that the BDD model 
mimics the wetting model in all respects that may affect universal scaling behavior. 

A clear insight into interfacial scaling as one approaches the scaling regime is obtained 
from Fig. ^(i) which shows the saturation width against system size for a range of desorption 
parameters k. With both axes logarithmic, one can define an effective exponent a{L, k) 
as the gradient of the function \ogW{\ogL, k). As is evident from the graph, there is a 
significant range of L, k for which a(L, k) 0.4-0.45 (similar to the values of P]). However 
the graph also gives evidence for a trend in which the limit lim^^oo (^{L, k) coincides for all 
values of k. This limit would give the true scaling value for the roughness exponent a for 
all K. We show this by collapsing all the data of Fig. ^(i) onto a single curve. 

To effect the collapse we use the simple crossover scaling form 

W,^tiL,^) = L''[A + Bie{^)/Ly] (4) 

where £{k,) is some finite length induced by the nonlocal desorption process and 7 is the 
crossover exponent. Remarkably, a reasonable collapse for all L, k could be effected by 
taking A,B constant and £{k,) = expXn with A constant — see Fig. ^(ii). That is, all the 
K-dependence enters through £(k) in a very simple way. 

The fit to (P allows a precise estimate of a=0.50(l) (coincident with the KPZ value) 
for the BDD model for all n. Without invoking (^ the estimation of a would be hampered 
by the slow power-law crossover to the asymptotic regime and underestimates of a would 
be obtained (as discussed above). This crossover would explain the discrepancies between 
previous estimates of a for BD (k=0) and the KPZ value P,§J^. 

It is interesting to see how the crossover scaling form (H) generalizes the standard pro- 



cedure of introducing an intrinsic width, first used for an Eden growth model [|T^,|T^. In 
that approach one writes W^^^ = AL'^°' + wf, i.e. the scaling width and the intrinsic width 
Wi are added in quadrature. Our scaling form would coincide with this procedure for large 
L if 7 = 2q; = 1. Moreover, 7 = a = | would correspond to the distinct procedure of the 
linear addition of two widths. For the BDD model, we have found that data collapses, all 
of reasonable quality, could be achieved for values of 7 ~ 0.5-0.7. 

To summarize, we first studied a wetting model whose bulk phase contains holes if a 
parameter k is nonzero. In the case k=0, the expected KPZ scaling was observed whereas 
for k nonzero simulations on small systems gave a value of the roughness exponent a rem- 
iniscent of values reported for a wide range of models with holes in the bulk phase P,^. 
Through theoretical considerations we concluded, in fact, that with k nonzero and in the 
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FIG. 4. (i) Log-log plot of saturation width Wgg^^ for different L and n in the BDD model. 

The solid lines correspond to a=0.5 (KPZ), 0.4 and 0.3 and illustrate how underestimates in a 

arise if one has not taken care to reach the scaling regime, (ii) The same data scaled according to 

Eq. (0). The parameters a=0.5 and A=7.4 were used and the solid line corresponds to ^=0.47, 

5=0.66 and 7=0.55 in @. 
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true asymptotic scaling regime, KPZ exponents should return. To demonstrate and quan- 
tify the crossover to the KPZ regime we studied the BDD model for which more extensive 
data could be obtained. We identified a power-law crossover of the form (^. An interesting 
question is as to the applicability of this simple form to other models, and in particular 
whether the value of 7 is universal. 

To conclude we place our observations in the context of other work. In a very recent 
preprint []TU| it has been suggested that in the case of a stable phase invading an unstable 
phase (e.g. Fisher waves and our wetting model) one should observe not the d+1 KPZ 
exponents but instead those of the KPZ equation in d+2 dimensions. 

The essence of that work is that in the situation where an interface is not sharp (e.g. due 
to the presence of holes), the correct surface to consider in terms of the scaling is (a suitable 
transformation of) the density field interpreted as a height variable. As the dimensionality 
of the density field over the smeared-out interfacial region is necessarily one greater than the 
interface, the d dimensional interface should scale in the same way as a line across the d+1 
dimensional density surface governed by the d+2 KPZ equation. However, as also pointed 
out in |T0[, the system does not scale isotropically: a rescaling transformation would affect 
only the size of the substrate. Thus if the interfacial region remains finite for increasing 
L the relative size of the extra dimension shrinks to zero and thus one returns to the d+1 
KPZ equation, i.e. there would be a crossover to KPZ scaling. For the models studied in the 
present work, our numerical evidence explicitly shows a finite interfacial region indicated by 
the typical gap size remaining constant as the system size L is increased. Furthermore, for 
the BDD model we could quantify the crossover. It would be interesting to confirm whether 
this crossover is also present in the model of p]. 

We thank Alastair Bruce and David Mukamel for helpful suggestions and EPSRC for 
financial support (RAB). 
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